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ABSTRACT 

We consider the physics of free precession of a rotating neutron star with an 
oblique magnetic field. We show that if the magnetic stresses are large enough, 
then there is no possibility of steady rotation, and precession is inevitable. Even 
if the magnetic stresses are not strong enough to prevent steady rotation, we 
show that the minimum energy state is one in which the star precesses. Since the 
moment of inertia tensor is inherently triaxial in a magnetic star, the precession 
is periodic but not sinusoidal in time, in agreement with observations of PSR 
1828-11. However, the problem we consider is not just precession of a triaxial 
body. If magnetic stresses dominate, the amplitude of the precession is not set 
just by the angle between the rotational angular velocity and any principal axis, 
which allows it to be small without suppressing oscillations of timing residuals 
at harmonics of the precession frequency. We argue that magnetic distortions 
can lead to oscillations of timing residuals of the amplitude, period, and relative 
strength of harmonics observed in PSR 1828-11 if magnetic stresses in its core 
are about 200 times larger than the classical Maxwell value for its dipole field, 
and the stellar distortion induced by these enhanced magnetic stresses is about 
100-1000 times larger than the deformation of the neutron star's crust. Magnetic 
stresses this large can arise if the core is a Type II superconductor, or from 
toroidal fields ~ 10 14 G if the core is a normal conductor. The observations of 
PSR 1828-11 appear to require that the neutron star is slightly prolate. 

1. Introduction 

The convincing observation of free precession of PSR 1828-11 (Stairs, Lyne & Shemar 
2000) poses challenges for theories of neutron stars. Shaham (1977; 1986) argued that 
vortex line pinning in the neutron star crust should prevent long term precession. Sedrakian, 
Wasserman & Cordes (1999) showed that precession is still prevented if vortex lines are 
not pinned perfectly but vortex drag is strong. They also showed that even if vortex 
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drag is weak, precession is damped away. Link & Cutler (2002) estimated the strength 
of vortex line pinning forces, and argued that PSR 1828-11 may precess at large enough 
amplitude to unpin superfluid vortices in the crust. If vortex drag is small this would 
remove one impediment to free precession, although the free precession would still damp 
away eventually if it is not excited continuously. 

Here, we consider an additional feature of radiopulsars like PSR 1828-11, namely, that 
they are strongly magnetized, with magnetic axes that are at an angle to their rotation 
axes. Based on earlier work on magnetic stars by Mestel and collaborators (Mestel & 
Takhar 1972, Mestel et al. 1981, Nittman & Wood 1981; see also Spitzer 1958), we argue 
that precession may be required - there is no equilibrium corresponding to solid body 
rotation without precession for a rotating star with an oblique magnetic field. For a fluid 
star, though, we shall see that although the fluid must precess, the magnetic axis rotates 
uniformly. Although Mestel et al. (Mestel & Takhar 1972, Mestel et al. 1981, Nittman 
& Wood 1981) showed that hydrostatic balance also requires fluid motions in addition to 
the precession which affect the stellar magnetic field, these are too slow to be important 
observationally. 

Once we also take account of the solid crust of a neutron star in addition to its fluid 
interior, we show that if the stellar distortions due to the magnetic field are larger than 
the distortion of the crust, then steady state rotation is very unlikely (but not necessarily 
impossible). If the magnetic stresses inside PSR 1828-11 are simply due to the classical 
Maxwell stress tensor, evaluated with the inferred dipole magnetic field strength, then they 
are too weak to require precession at a period ~ 1 year. However, if the interior of the 
neutron star consists of a Type II superconductor, the effective stress tensor is larger than 
for a classical magnetic field, according to Jones (1975) and Easson & Pethick (1977); as was 
emphasized by these authors, and Cutler (2002), the magnetic distortion is correspondingly 
larger. For PSR 1828-11, we estimate that the distortion that would result in a core that 
contains a Type II superconductor can lead to precession at a period of order 1 year. A 
sufficiently strong toroidal magnetic field (B t ~ 10 14 G) could also lead to precession at this 
period without Type II superconductivity (see e.g. Eq. [2.4] in Cutler [2002] with B c — > 0). 

Even if magnetic stresses are not strong enough to prevent steady rotation, magnetic 
distortion in an oblique rotator alters the physics of free precession qualitatively and 
quantitatively compared to what one would expect for precession due to axisymmetric 
crustal distortions. Even if the crust is axisymmetric, misalignment between the symmetry 
axis of the crust and the magnetic field make the effective stellar moment of inertia 
inherently triaxial. From a phenomenological viewpoint, one manifestation of this loss 
of axisymmetry is that although the angular velocity of the star is a periodic function 
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of time in the frame rotating with the crust and magnetic field, it is not a sinusoidal 
function of time. This feature is consistent with observations of PSR 1828-11, which reveal 
behavior at several different harmonically related frequencies (Stairs, Lyne & Shemar 2000). 
The relative strengths of the harmonics depend on the degree of nonaxisymmetry, and, 
presuming an axisymmetric crust, on the distortions induced by the magnetic field. For 
the ordinary Maxwell stresses evaluated just with the inferred dipole field strength, the 
nonaxisymmetry would be small, but distortions resulting from a Type II superconductor, 
or from a normal core with a large toroidal field, could produce sufficient nonaxisymmetry 
to lead to comparable amplitudes for at least the first few harmonics of the fundamental 
precession period, as is observed. 

Periodic, but not sinusoidal, precession would also arise if the neutron star crust were 
simply nonaxisymmetric even if magnetic stresses were negligible. Thus, the detection of 
harmonic behavior in PSR 1828-11 cannot, by itself, be taken to be evidence for amplified 
magnetic stresses in its core. However, standard results for triaxial precession, which are 
reproduced as a byproduct of the calculations we present, show that the oscillations at 
harmonics of the precession frequency are smaller in amplitude than the oscillation at the 
fundamental frequency by powers of the precession amplitude. We shall see that this is not 
the case in models where magnetic stresses predominate: the amplitudes of oscillations at 
the precessiom frequency and twice the precession frequency may be comparable even at 
small amplitude. 

Moreover, it is well known that the minimum energy state for rotation of a 
nonaxisymmetric body is one in which the angular velocity and angular momentum are 
aligned with the principal axis of the body with largest eigenvalue of the moment of inertia 
tensor. There is no precession at all in this state. However, we shall see that even when 
magnetic stresses are not strong enough to require free precession, the minimum energy 
state for an oblique rotator does not correspond to alignment of the angular velocity with 
any principal axis of the effective moment of inertia tensor. Thus, the minimum energy 
state is one in which the star precesses. 

We shall see that the timing residuals associated with the precession can account 
for observations of long term oscillations in PSR 1828-11, but that the explanation only 
works if the core of the neutron star has sufficiently strong magnetic stresses that magnetic 
distortions are 100-1000 times larger than crustal deformations. Thus, we argue that 
observations of free precession in PSR 1828-11 offer evidence that the pulsar is in the regime 
where magnetic stresses dominate, because either its core is a Type II superconductor, or, 
if it is a normal conductor, possesses a very large toroidal field. 

We review the arguments given by Mestel and collaborators (Mestel & Takhar 1972, 
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Mestel et al. 1981, Nittman & Wood 1981) in In §[| we consider the conditions 
under which precession of an oblique rotator is inevitable. In §[|, we consider the modified 
precession problem for an axisymmetric crust and misaligned magnetic field. There we 
show that the minimum energy state at a given angular momentum is one in which the 
star precesses, provided that the magnetic field is not either along or perpendicular to the 
symmetry axis of the crust. There, we also review classical results for free precession of 
triaxial bodies; we shall find that there are three distinct cases of interest for a star where 
magnetic distortions are important. We present limiting results for the timing residuals 
expected in this model in §[|, and obtain approximate results for the limit in which crustal 



distortions dominate in §5.1, and for the opposite limit in which magnetic distortions 



dominate in 85.2. 



Here, our main purpose is to present arguments that an oblique rotator must precess. 



In § |5.3| , though, we present brief application of the model to observations of PSR 1828-11. 
There we argue that only models in which magnetic distortions dominate can account for 
all of the observed features of the long-term periodic timing residuals from this pulsar. We 
also suggest that the data favor prolate rather than oblate magnetic distortions. 

Several appendices present cumbersome mathematical details needed to derive (and 
verify!) analytic results presented in the text. 



2. Review of Structure of Fluid Equilibria with Oblique Magnetic Fields 

Let us begin by reviewing the theory of rotating stars with oblique magnetic fields 
developed by Mestel and collaborators (Mestel & Takhar 1972, Mestel et al. 1981, Nittman 
& Wood 1981; see also Spitzer 1958). Consider a (fluid) star with angular momentun 
L = LL. In the absence of a magnetic field, the angular velocity of the star would be 
f2o = L/I = Q Q L. Let us work in a reference system rotating at a rate u> = uuL, where u is 
unknown, but must be Qq to lowest order in small quantities. Assume that the magnetic 
field is axisymmetric about an axis b that is fixed in this frame; we can verify that this 
works out correctly to lowest order later. 

Solve for the density field of the star under the assumption that it is stationary; Mestel 
et al. (Mestel & Takhar 1972, Mestel et al. 1981, Nittman & Wood 1981) discuss the sizes 
of time dependent correction terms, which are smaller than any we retain here. The result 
is a superposition of three components assuming slow rotation and weak magnetic fields: 

p(r) = p (r) + p u (r, r-L) + p B (r, f-b) . (1) 

Here, po( r ) is the spherical profile of the undistorted star, and p w (r, f-L) and ps(r, f-b) 
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are the distortions due to rotation and magnetic fields, axisymmetric about L and b, 
respectively. In absolute magnitude, the rotational distortion is of order u 2 R 3 /GM for a 
star of radius R and mass M, and the magnetic distortion is of order BHR^/GM 2 , where 
H = B for ordinary magnetic fields, but H corresponds to the first critical field strength in 
a Type II superconductor (H c \ ~ 10 15 G), as discussed by Jones (1975), Easson & Pethick 
(1977) and Cutler (2002). 

The moment of inertia tensor corresponding to Eq. ([]]) is of the form 

Iij = J d 3 r p(r) (r 2 5ij - r^) = J 5 y + I u - 3£rjj[y) + I B (<% - 36*6^) , (2) 
where (P 2 (/i) = \ (3p 2 - 1)) 

2 /" 

h = -J d 3 r p(r)r 2 

L = \j 'd 3 r p„(r,r-L)P 2 (r-L) 

Ib = \j d 3 r p B (r,r-b)P 2 (r-b) . (3) 
Let ft = L/I + 8Q = Qq + 5f2; then the angular momentum of the star is 

L = m = {I Q +I w +I B )n-MuLL-n-M B bb-n ~ (I -2J w +J B )n +I<$n-3I B b&-no , (4) 
which implies that 

I y t~tJ — 1 — ■ 

We can decompose this into components along and perpedicular to L i.e. <5f2 = 5Q,\\L + 5Q±, 
where 



l-^ + ^ + ^P 2 (b.L) 
III ' 



<JJ2|I ~ fi c 

<5fi_i_ = only if b is either parallel to or perpendicular to L. Since Euler's equations imply 
time independent Q only if L and f2 are parallel, in general a rotating star with an oblique 
magnetic field must precess. 

To find the angular velocity of the reference system, require that the magnetic axis 
rotates with angular velocity fi = uiL + 3bI B b-fl / P, thus 

« = n.(2-£ + ^-£) (7) 
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In this reference system, matter precesses about b with angular velocity 

3i B n bb-t 

Up = J • (8) 

However, even though the fluid precesses, the magnetic field does not, so we do not expect 
any observable effects to arise from the precession. 



3. Conditions Under Which Steady Rotation is Impossible and Precession is 

Inevitable 

We can look at the Mestel et al. (Mestel & Takhar 1972, Mestel et al. 1981, Nittman 
& Wood 1981) results in a slightly different way, which is identical to the same order of 
approximation, but differs in approach slightly. Work in a reference system corotating with 
the matter. [] The perfect conductivity condition demands that the magnetic axis is fixed 
in this reference system. Instead of Eq. (H) , suppose that 

p(r) = p (r) + p n (r, f-Cl) + p B (r, f-b) ; (9) 

in this case, the moment of inertia tensor (formerly Eq. J|) becomes 

kj = J d 3 r p(r) (r% - nr 3 ) = J c% + I a fa - 30,0,) + I B fa - 36^) , (10) 

where the various factors are defined just as in Eq. (|3|). The angular momentum is therefore 

L = (I -2I n + I B )n-3I B bb-n. (11) 

The condition for a time independent Q is L \\ O. This condition is satisfied only if b is 
either parallel to or perpendicular to Q ~ L. The star precesses at an angular frequency 
uj D p as before. However, the magnetic axis rotates at a uniform angular velocity, because, 
from Eq. ([□]) 

Slxb= L * h =Q B xb, (12) 
Jo — + Ib 

where 

n B = ~ o (2 - ^ + ^ - cj , (13) 

I Q + 2I a + I B °\ I I I J K } 



Mestel et al. (Mestel & Takhar 1972, Mestel et al. 1981, Nittman & Wood 1981) show that maintaining 
hydrostatic balance also requires additional motions of the fluid. We shall discuss this feature in more detail 
below. 
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which is independent of time. 

The second approach makes contact with the Euler problem clearer. Angular 
momentum conservation is consistent with time independent rotation as long as the angular 
velocity vector is aligned with one of the principal axes of the moment of inertia tensor. For 
anisotropic density distribution, the moment of inertia tensor is of the form 

hj = lo^ij + A/jj 

where Tr(Aiy) = 0. Even though ||Wy|| <C Iq for small distortions, the principal axes of /y 
are the principal axes of AI^. For the magnetic fluid, Eqs. (||D and (0) imply that 



Alij = I n [6 tj - ZSliSlj) + I B [Sij - 36 i 6 j 
For there to be a principal axis of Ay along f2, 

= AO, . 

where A is the associated eigenvalue. For this form of Ay, the condition becomes 

-3I B bib ■ O = (A + 2/ n - I B )Sl u 

which is only true in general provided that either ft is along b or O is perpendicular to 6, 
neither of which will be the case generally. 

Next, consider what happens when we consider a neutron star model that consists of a 
rigid solid and a fluid, with an oblique magnetic field. In that case, Eq. fllOD is generalized 
to 

hj = lot* + + 1^ - 30^-) + I B (Sij - Sbibj) , (14) 

where is the trace-free part of the moment of inertia of the crust. The angular 
momentum of the star is therefore 



Li = (Jo - 2/ n ) ft, + [Cij + I B [Sij - 3^6,] j Sl d . (15) 
A steady state is only possible if the star rotates along a principal axis of 

Alij = C lj + I B [S^ - 3bibj 
But for this to be true, we must have 

- 3bj B {b-n) = (A - i B )n t , (16) 



2 More precisely, Cij is the portion of the moment of inertia tensor of the crust that is not aligned with 
the magnetic distortions. Any crustal distortions symmetric about b would just renormalize 1b- 
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where A is the associated eigenvalue. The components of this equation perpendicular to O 

are 



Eq. (|T7D must have a solution (with Cl — L) in order for there to be no precession. 



;i7) 



To understand the significance of Eq. (|17|), first review the argument against precession 
when magnetic fields are ignored. The moment of inertia tensor of the solid is not known 
observationally for any neutron star, so we can specify freely. Since it is a trace-free 
symmetric tensor, Cy has five independent components. We only need to choose two of 
these (the orientations of one of the three orthogonal principal axes) to get a steady state. 
We can make these choices without considering the magnitude of the distortion of the solid, 
since only the directions of the eigenvectors matters for finding a steady state. Moreover, 
from a physical viewpoint, we expect that at a fixed L the energy associated with rotation 
is minimized if L (and hence SI) is along the principal axis with the largest eigenvalue; 
dissipative processes tend to drive the neutron star toward this state. 

Next, let us consider what happens when we restore the magnetic field. Expand 



— C^Aj A 



AC/ 1 

i 



where are the eigenvalues and A M the eigenvectors of Cy, and rewrite Eq. ([T7|) as 



a^a" 



3I B (b-n) 



Q(b-Q) 



(19) 



We know that if I B = 0, a steady state is possible irrespective of the magnitudes C„, and 



that if all C M = there is no steady state. When I B is nonzero, we can adjust the A M so 
that Eq. fll9"|) is satisfied, just as for I B = 0, but the adjustment depends on the magnitudes 
C M relative to I B . In general, we would expect there to be no solution if |C^| ^ Ib, in which 
case the neutron star must precess. In § |4] we shall also see that even when I B is not large 
compared to C, the rotational energy is minimized, to first order in distortions, when the 
angular velocity is not precisely along a principal axis of the effective moment of inertia 
tensor, and therefore not precisely along L. 

To appreciate this result better quantitatively, let us focus on the simplest special case 
in which the crust is axisymmetric. If k is the symmetry axis of the crust, then we may 
write 



a 



(20) 



so that the largest eigenvalue of the moment of inertia of the crust is along k. Then the 
condition for a steady state, Eq. (|19|), may be written as 



3C(k-{l) 



k - rt(k-n) 



3I R (b-Q) 



n(b-n) 



o . 



(21) 
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If Eq. ( pip has a solution, it must have coplanar k, bCl. Thus, let us adopt b = e s and 

fl = e 3 cosx + ei sinx 

k = e 3 cos 9 + ei sin 9 : (22) 

substitute into Eq. (^l|) to find 

(e 3 sin x — &i cos x) [C cos(x — 9) sin(x — 9) — Ib sin x cos x] = . (23) 
Therefore, for a steady state, we must require 

S in[2(x- e) ] = ^, (24) 

which is only possible if |/ssin2x| < C. If Ib ^ C, then steady state rotation is rather 
unlikely, with a probability that decreases with increasing Ib/C. When Ib sin2x/C < 1, 
Eq. (plf) has the solutions 



/Bsin^x , _ / /lsin 2 2 X 
cos 29 = ±cos2xyl — 2 



. . , / J|sin 2 2x Jbcos2x\ . . 

sin 20 = sm2x|±yl ^ ^ , (25) 

as can be verified by direct substitution; the choice of signs depends on specific parameter 
values. 

In order of magnitude, 3/^/1 = (3BHR 4 / 'GM 2 , where (3 ~ 1 is a structure constant, so 

that 

lob = 3IbQo i COSX = 1.9 x 10-"n (3 cos X (BH) 27 RtM^ , (26) 

for a neutron star mass and radius M = 1.4M 14 M and R = IORq km, a magnetic axis 
inclination angle x = cos" 1 (&•£), and BE = 10 27 \BH) 27 G. (Si? ~ 10 27 G for B ~ 10 12 G 
and if ~ 10 15 G ~ i7 c i - see Easson & Pethick 1977.) The period associated with ojb is 
Pb = 2h/uj b — 16 .6 Po(s) (f3 cos x)~ 1 M 2 a Rq 4 (BH) 2 7 y for a neutron star rotation period 
P = lPo( s ) s - If magnetic effects are strong enough to require precession, then we should 
expect a precession period of order Pb- We consider the precession period more completely 
in the next section. 

For PSR 1828-11, the spin period is Pq = 0.405 s and the dipole field strength deduced 
from the observed spindown is B ~ 5 x 10 12 G (Stairs, Lyne & Shemar 2000); these values 
would imply P B — 1.35(/?cosx) _1 [5/(j9ii')27] y, or about 492(/3cosx)~ 1 (-Bif 2 7/5) d, similar 
to the observed period for flcosx ~ 1- Note that the precession period would be far longer 
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for BH — > B 2 , by a factor of about 200. Thus, if magnetic effects are the reason for the 
observed "precession" then ether the neutron star interior must be a Type II superconductor 
(Jones 1975, Easson & Pethick 1977, Cutler 2002), or there must be a substantial toroidal 
field in the core if it is not a Type II superconductor (e.g. Eq. [2.4] in Cutler 2002 with 
B c ^0). 



4. Modified Euler Problem 

4.1. Basic Equations, Principal Axes and Eigenvalues 

Next, let us consider the modified Euler problem. To allow an analytic treatment, 
continue to assume that the crust is axisymmetric. Although this is a simplification, 
the effective moment of inertia of the star will still be triaxial because of the misaligned 
distortion introduced by the magnetic field, except for special orientations. 

The Euler equation in the rotating frame of reference is 

/;/ 7 ^ + 3ey,fy (Ckk m - I B bib m ) On = , (27) 
where the effective moment of inertia tensor is 

if = (J - 2I n + I B - C)5 tJ + 3 (Ckikj - I B bib } ) . (28) 

and d*/dt is the time derivative in the rotating frame. The two vectors, b and k define a 
plane, and one of the eigenvectors of is along the unit vector perpendicular to that 
plane, e 2 , and has an eigenvalue J 2 = J — 2Iq + I B — C. The other two eigenvectors lie in 
the b — k plane. As in the previos section, take £3 = b and define k by Eq. (|22|). Then the 
other two eigenvectors, e± are 3 

e± = e 3 cos <r± + ex sin a± , (29) 

where 

1 ( v V /2 

sm a± = —= 1 ± . 
3 The following results imply that 

e+ xe_ = &2 
e 2 Xe + = e_ 
— e 2 Xe_ = e + 
so that these axes define a right handed coordinate system. 
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cos a± 



V 



I B -C cos 26 



1/2 



Csin2# 

The eigenvalues associated with these eigenvectors are I± = I2 + where 



C-I B ± \ C 2 



2CI B cos 29 + II 



(30) 



(31) 



In Appendix |A| we show that for Is > 0, i+ > I2 > I-, and for Ib < 0, /+>/_> I 2 . 

Approximate results are derived for \Ib\ <C C in Appendix [B| and for \Ib\ 3> C in 
Appendix |C]. When \Ib\ C, the eigenvectors are nearly aligned with the principal axes 
of Cjj, and the eigenvalues are nearly the eigenvalues of CV,-, apart from small corrections 
~ \Ib/C\. For \Ib\ 3> C, the eigenvectors are nearly along and perpendicular to b and the 
eigenvalues are almost determined by the magnetic distortions alone, apart from corrections 
~ \C/Ib\- We shall argue below that the case Ib < 0, \Ig/C\ ^> 1 may be especially 
relevant to observations of PSR 1828-11. 



4.2. Solution of the Euler Equations 

4-2.1. Basic Equations and Minimum Energy State 
The modified Euler equations are simply what one finds in general for a triaxial system, 



si-SizSi- = 



^ alt 

i^—^ + 5i + n 2 n+ = 

at 

- ( 5I + ~ = . (32) 

As is well-known, Eqs. ([32]) conserve both the magnitude of L and the rotational energy, 

E rot = I (i + n 2 + + i^n 2 _ + i 2 n 2 2 ) . (33) 



Note that the magnitude of the angular velocity is not conserved, so that in actuality Iq is 
not independent of time. However we will assume that the variation is slow enough, and 
of small enough amplitude, that its effect is only higher order than any others we consider 
here. 
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The condition for steady rotation, Eq. (|24|), can be rederived from Eqs. ([[(]) under the 
assumption that Cl is along either e + or e_, with \ < 7r /2 or % > 7r/2, respectively. We 
are interested in what happens when precession is required, so let us assume that Eq. (^|) 
cannot be satisfied for any choice of 9 given x, Ib and C. Nevertheless, we expect that the 
system seeks a minimum energy state, even if that state is not steady. To be specific, let 
X be the angle between b and Cl when they are coplanar with k (and therefore with e±). 
Let us determine 9 from the requirement that E is minimum for a given L 2 . Thus, let us 
consider the quantity 2E/L 2 at the epoch when Cl, k and 6 are coplanar. If we define 



/ = 
51 = 

A = 



(/+ + /_)=/ 2 + 



3(C-/ B ) 



2 
51 

T 



(A 



2CI B cos29 + I 2 B 



2I 2 + 3(C - 7 f 



then we find 



2J£ 



i + A(2ri + - 1) 



L2 1 + A 2 + 2A(2f2 + - 1) 
where the approximation holds for A< 1, and 



+ 



V2~{ 



cosx 



1 



A(2f2+ 



(I B -Ccos29) 



+ sinx 



1 + 



C 2 

(I B -Ccos29) 



1/2 



2CI B cos29 + I B 

1/2 



ft 



smx 



2CI B cos29 + I 2 B 
(I B -Ccos29) 



1/2 



cosx 



2CI B cos29 + C 2 

1/2 



1 + 



(I B -Ccos29) 



Assuming that 9 < n/2 we find 



A(2Ct_ 



3{Ccos[2( X 



2CI B cos29 + C 2 



I B cos 2x} 



2I 2 + 3(C - Ib) 



(34) 



(35) 



(36) 



(37) 



To first order in A, the energy is minimized when A(2f2 
when 9 = X- 



1) is maximized, which happens 
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Note that in the minimum energy configuration, 0_ need not be very small, although 
it vanishes in the minimum energy state for I B — > 0. For small values of Ib/C, 

Cl ~ sin( X sin 29 cos( X - 0) -> , (38) 

where the last result assumes 9 — > x- Thus, for small is, we find small but nonzero Cl- 
in the minimum energy state. For \Ib/C cos2#| ^> 1, Cl ~ — sin 2x/2C. On the other 
hand, for \I B /C cos29\ 1, Cl ~ sin% or — cosx, depending on the sign of Ib- 

When I B = 0, 9 = x minimizes 2EI/L 2 exactly, and we find, as usual, + = 1, 
corresponding to an angular velocity aligned with the principal axis with the largest 
moment of inertia. When I B ^ 0, 9 = x does not correspond to exact alignment of the 
angular velocity and the principal axis with the largest moment of inertia. This is because 
the eigenvalues of 1^ depend on 9, so A depends on 9. If A were independent of 9, then 
the energy would be minimized for Cl + = 1, the largest possible value of Cl + . The value of 

LCl = - ! + K ! _ 2A 2 Cl 2 + Cl 2 (39) 

VI + 2A(2^ + - 1) + A 2 

is only exactly one when either Cl + = 1 or 0_ = 1; the angle between L and Cl is 
±2Af2 + f2_ more generally. 

The condition that 9 = x means that the angle between the angular velocity vector 
and the magnetic field is the same as the angle between the symmetry axis of the crust and 
the field. Thus, it corresponds to angular velocity that is parallel to the symmetry axis 
of the crust at the epoch when all of the vectors lie in a plane. Since the angle between 
the angular velocity vector and the magnetic axis is fixed in this case, the portion of the 
rotational energy due to the rotating magnetic distortion is fixed. Minimizing the energy 
then amounts to minimizing the portion of the energy associated with rotation of the 
crustal distortions. This is achieved if the angular velocity is along the symmetry axis of 
the crust. The resulting minimum is only a local minimum, at a given value of x- The 
global minimum is achieved for x — or ir/2. We assume that the star can evolve on a 
slow disspative timescale toward the local minimum, and on a longer timescale, probably 
the spindown timescale (e.g. Goldreich 1970) toward the global minimum. 

The fact that A = A(cos2#) when Ib ^ implies that even when alignment of O 
and L is possible, the minimum energy state is not steady rotation, but rather precession. 
When I B C, the minimum energy state corresponds to an angle of approximately 

~ \Ib sin2x|/2C between Cl and e+, and an angle ~ 3\Ib sin2x|/2J between Cl and L. 
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4-2.2. Nonlinear Solution: Ib > 



The complete, nonlinear solution to Eqs. ( p2[ ) is given in Landau & Lifshitz, §37, in 
terms of elliptic functions. We will have to consider separately the two cases Ib > and 
Ib < 0. Here, we consider Ib > 0, which implies I + > I2 > I according to Eq. (\K%) and 
the ensuing discussion. Adapting the solution in Landau & Lifshitz, §37, to this situation 
(and our notation) we have (see their Eqs. [37.8]- [37. 12]) [] 



L 



fi_cn(Y) 



1 + 2A(20 + - 1) + A 2 
f2_sn(r j 



2A(1 - A) 



\ [1 + 2A(2f2 + - 1) + A 2 ](l - A )(A + A c 
f2+ dn(r) 



1 + 2A(20 + - 1) + A 2 



Q+dn(r) 



(40) 



where A = 3(C-I B )/2I, cn(r) = Jl - sn 2 (r), dn(r 



sii(t) 



dx 



with 



yj \l ~ X 2 )(l - k 2 X 2 ) 



(A - a )(i - A)(i - n+) (A - a )(i - n 2 + 



1 — /c 2 sn 2 (r), sn(r) is defined by 

(41) 



-2 



dr 
~dt 



L 



(A + A )(l + A)ri 



(A + A )^+ 



2A(1 + A)(A + A 



1 \ (1- A )(l- A 2 )[l + 2A(2f2 + - 1) + A 2 



The motion is periodic, with a dimensionless period 4K(k 2 ) where 

dx 



K(k 2 



J(l-x 2 )(l-k 2 x 2 ) 



(42) 



(43) 



4 Thcrc is actually a sign ambiguity in the solution given by Landau & Lifshitz, which we resolve by always 
choosing f2_ oc fl , rather than £1_ oc *J 1 — fl + . The Euler equations Eq. ( p7| ) require that if wc choose 
57_ oc f2_, then we should also choose the sign of the coefficient of Vl 2 to be the same as the sign of the 
coefficient of 
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Eqs. (f40|) and (|42]) also imply that 



I 2 W _ 1 - [2A(A - A )(l - ft 2 + )/(l + A)(l - Ap)] sn 2 (r) 
L 1 + 2A(2fJ + - 1) + A 2 

which is independent of time up to terms ~ A 2 , thus validating the approximation of time 
independent In used above. 

One of the distinguishing features of the timing model we will develop below is that we 
shall not demand that |f2_| be small. In particular, we shall see that when \I B \ 3> C, Cl_ 
will not be small in general, but the observable effect of the precession on pulse arrival times 
could still be small. This is because in the limit where magnetic distortions are far larger 
than crustal distortions the star tends to precess about its magnetic axis. If there were no 
crust at all, as in the magnetic fluid stars considered by Mestel and collaborators (Mestel 
& Takhar 1972, Mestel et al. 1981, Nittman & Wood 1981), the star would precess exactly 
around its magnetic axis, which would therefore rotate uniformly. The crust breaks this 
symmetry, and allows the precession to be observable. We shall see this emerge in some 
detail when we consider timing residuals in §|5~2\ 



For the opposite case, where the crustal deformations dominate, the precession 
amplitude is set by 0_, according to Eqs. (|40|). In that case, we also see that Eq. (|42|) 

~ 2 

shows that k 2 oc f2„. Thus, if the precession amplitude is small, so is k 2 irrespective of how 
triaxial the crust might be. Since k 2 governs the importance of oscillations at harmonics of 
the fundamental precession frequency, we see that small amplitude will imply oscillations 
predominantly at the fundamental if crustal deformations dominate. (This will also be true 
for the solution given in §[4.2.3| for I B < 0, where Eq. will also imply that oscillations at 



harmonics of the precession frequency are suppressed for small precession amplitude.) This 
is not the case when magnetic deformations dominate. 

The better-known results for free precession of an axisymmtric star are recovered for 

A — A = = k 2 . However, note that in the triaxial case, k 2 need not be small. We have 

- 2 

already noted that 1 — f2 + does not have to be small, even in the minimum energy state, 
and from Eqs. (B3) and the definition of Aq we see that 



A-A JC 2 -2CI B cos26 + I 2 - {C - I B ) , * 

= — (45) 

A + A o ^C 2 - 2CI B cos 26 + I 2 + C - I B 

which is not necessarily small either. For < I B /C <C 1, we note that 

k 2 ~ I B sin 2 2% sin 2 x/4C 3 to lowest order in I B /C in the minimum energy state, so 

k 2 1 in this case. However, for I B ^> C, we shall see that k 2 ~ I B /C ^> 1, and in general, 
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since A < for Ig > C, it is possible for k 2 to exceed one rather generally in that regime. 
The solution to the Euler problem is still given in terms of elliptic equations when k 2 > 1, 
but we have to make the replacements 



sn(r) 
cn(r) 
dn(r) 



sn(fcr) 
k 

dn(fcr) 
cn(kr) 



(46) 



in Eqs. 



L 



where sign(fL 
and 



and dn 2 (/cr) = 1 — k 2 sn 2 (kr). The explicit solution for k 2 > 1 is 
n_dn(f) 



1 + 2A(2fi + - 1) + A 2 
+ sign(f2_)sn(f 



2A(1 + A) 



\ (1 - A )(A - A )[l + 2A(2f2 + - 1) + A 2 ] 



+ cn(r) 



— , (47) 
'1 + 2A(2^; - 1) + A 2 

is the sign of f2_, which may be positive or negative, dn 2 (f ) = 1 — /c 2 sn 2 (f), 



dr 
~dt 



2A(A - A )(l - n. 



I \ (1 - A )(l + A)[l + 2A(2f2 + - 1) + A 2 



k 2 = -r~ 



1 



(A + A )(l + A)S7 H 



(48) 



(A-Ao)(l-A)(l-n + ) 
This somewhat unfamiliar case is not treated in Landau & Lifshitz, but is found eaasily 



from the equations given there, and corresponds simply to the transformations in Eqs. (|4"6j ). 
Physically, it turns out to be important for a substantially prolate figure, which is what 
happens when Ib is large and positive. 

An important difference between axisymmetric and triaxial precession is that even 
though both are periodic, only the axisymmetric case is precisely sinusoidal. The functions 
cn(r) and sn(r) contain all odd harmonics of irT/2K(k 2 ), whereas the function dn(r) 
contains all even harmonics. When k 2 is small, the expansions are dominated by their 
leading terms, but, as we have seen, the general case does not demand small values of k 2 . 
The key parameter in the expansions is 



q(k 2 



exp 



TlK'(k 2 )' 

' K{k 2 ) 



(49) 



-17- 



where 

K\k 2 ) = K(l - k 2 ) (50) 

(Abramowitz & Stegun, Eqs. [16.1.1], [16.23.1-3]). Since q{k 2 ) < 1 by its definition, only 
the first few harmonics ought to be prominent in the solution, but as long as k 2 is not 
especially small compared to one, the amplitudes of the first few harmonics ought to be 
roughly comparable. This feature is consistent with the observed properties of PSR 1828-11 
(Stairs, Lyne & Shemar 2000). If this interpretation of the observations is correct, then the 
fundamental precession period must be 1000 (or perhaps 2000) days. f] 



4-2.3. Nonlinear Solution: Ib < 



For Ib < 0, the solution is analogous to what was given in §4.2.2j , except that, according 
to Eq. (|A3j) and ensuing discussion, I + > I_ > I 2 . Instead of Eqs. P0] ) we have 



L 

m_ 
m + 

and instead of Eq. 



ft 



2A(1 - A) 



1/2 



[1 + 2A(2ft + - 1) + A 2 ](l - A )(A + A ) 
f2_sn(r) 



cn it) 



1 + 2A(2fi + - 1) + A 2 



1/2 



dr 
~dt 



[A + A 2 + A(l + A )(2r2 + - 1)] 
[1 + 2A(20 2 + - 1) + A 2 ](l + A)(A + A ) 
j) we have 

(A -A)(l-A)(l-n+) 
A + A 2 + A(l + A )(2^-1) 

L 2A[A + A 2 + A(l + A )(2^ - 1)] 
7 \ (1-A )(1- A 2 )[l + 2A(2f2+-l) + A 2 



dn(r) 



(51) 



(52) 



The definitions of sn(r), cn(T) and dn(T) are the same as before. The magnitude of the 
angular velocity is now 



f\n\ 



1 + [2A(A - A)(l - n + )/(l + A)(l - A )]cu 2 (t) 



2A(2fi 



A 2 



(53) 



The longer period appears in the Fourier analysis of the timing residuals, but not in those of the period, 
period derivative, or shape variations, given in Fig. 2 of Stairs, Lyne & Shemar 2000. In any event, the data 
shown there only span about 2000 days, so this feature may not be so well-established. 



-18- 



Eq. (^) implies that |f2| is independent of time up to terms ~ A 2 , just as we found from 
Eq. (f44|) for I B > 0. In this case, k 2 < 1, so these solutions suffice. 



4.3. Effect of the Spindown Torque 



The solutions given in §4.2.2j and [4.2.3| are for free precession. Radiopulsars spin down 
as a result of electromagnetic radiation, on a timescale t s d = P/2P that is long compared 
with the precession periods of interest: for PSR 1828-11, t s d — 1.1 x 10 5 years. A careful 
examination of the Euler equations with spindown included shows that the precession is 
described by the torque-free solutions up to small corrections (just as was found by Link & 
Epstein 2001). We also note here that the quantity 2EI/L 2 only changes on a still longer 
timescale, ~ t s d/A. However, sinusoidal variations of the spindown torque result from the 
precession, since the angle between b and O varies with time. The amplitude of these 
variations need not be small, and the associated timing residuals can dominate (Cordes 
1993). In fact, we shall see that they are dominant, just as was found by Link & Epstein 
(2001) for precession of an axisymmetric neutron star. 



For vacuum magnetic dipole radiation, 



^ = Kbx(bxn), 

dt 



(54) 



where K is a constant to sufficient accuracy, and the magnitude of the angular velocity 
changes at a rate 



dfl _ldL 
dt I dt 



'it 



c 



sd 



sin 2 x. 



(55) 



where Qq is the angular frequency at some reference time t = 0; in the absence of precession, 
( = 0. We use the conventional definition of the spindown time as t s d = P/2P = —Q/2Q. 
We can use the solutions to the Euler problem found above to find the time dependent £ 
case by case: 



c 
c 
c 



O sn 2 M 



-2 



O , sn 2 (r" 



ft cn 2 M 



(A-Ao 


COS 2 <T + 


A + 


A 


(a + a ; 


sin 2 er + 


A- 


A 


(A-A 


) cos 2 er + 


A + 


A 



+ sin a + 



+ cos a. 



+ sin <7 + 



2f2 + fi_ cos<7 + sina + [l — dn(r)cn(r)] 
2f2 + f2_ sin a + cos a + [1 — dn(f )cn(f )] 



+2f2_ sin a + cos a + 



Ao + A(217 + -l) f 
dn (r sn r \ t-^ il 

y ' y >\ A + A 



(56) 
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for, respectively, I B > 0, k 2 < 1, I B > 0, k = 1/k 2 < 1, and I B < 0. Here, we do not 
restrict the solution to small values of f2„, and Eqs. ( [55]) and (|56| ) imply oscillations of 
the spindown torque at both even and odd harmonics of the precession frequency. (There 
are also evidently zero frequency corrections but these can always be combined with sin 2 % 
and factored out.) Note that for simple triaxial precession with a small tilt of the angular 
velocity away from the principal axis with maximum moment of inertia, the amplitude of 
the oscillations at 2uj p would be smaller, by a factor ~ \&-\, than the amplitude of the 
oscillations at u p . For the model developed here, f2_ can become substantial, particularly 
if I B < and \I B \ > C. 



5. Pulse Arrival Times 

To determine pulse arrival times, we need to represent b(t) in the inertial frame 
of reference of the observer. We can define the angular momentum vector L, which is 
conserved apart from spindown in this reference frame, to lie along the z axis, and we 
can further choose to place the observer in the x — z plane. Pulses arrive when b y (t) = 
(actually, only half of the solutions correspond to pulse arrival - the other half might be 
an interpulse or else unobserved) . The general problem is addressed in Appendix [D|, where 
the three different types of solutions are treated separately in Appendices [13. 1| and p.2| for 
positive and negative I B , respectively. Approximate solutions are also derived in those 
appendices (Eqs. |pi7|| , [p2^| and |D41| |), but those results only apply when the star is 



nearly axisymmetric and |f2_| 1. 

Here, we shall investigate two different nearly axisymmetric cases, \I B /C\ 1 and 
\Ib/C\ ^> 1. For \I B \ C, it will turn out that |f2_| 1, and the results of Appendices 
p.l| and p.2| will be directly applicable. However, for \I B \ ^> C, we have already mentioned 
that f2_ need not be small (see, for example, discussion following Eq. [BH]). 



5.1. Pulse Arrival Times for \I B \ C 

When \I B \ <C C, we can apply the results in Eqs. ( pi 7] ) and ( |D41| ) directly. Using the 
approximations given in Appendix [B], and accounting for pulsar spindown using the results 
of § [4.3|, we find that the oscillating parts of the phase residuals are 



- / Qq cos t \ 

il At osc ~ f2_ cot x - sin r + — — (57) 

V t sd ui 



p 
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for I B > 0, and 



for I B < 0, where 



il At osc ~ S2_cot % - cost + — 5- (58) 

V / 
n * - "+^2(9 (59) 

and fii'' = sin(% — 0), which is zero in the minimum energy configuration. (See Appendix 
|TJ.) Thus, to lowest order in the (presumed) small quantity f2_ cot x, phase residuals 
oscillate only at uo p . Oscillations at higher harmonics, such as 2u; p , have amplitudes that are 
smaller by additional factors of f2_. It is possible for these to be comparable in magnitude 



to the terms ~ f2_ retained in Eqs. ([57]) and (p8[) , but only if cotx is very small i.e. if 



- 2 

|0_ cot x| ~ or I cotx| ~ |f2_| <C 1. This was also found by Link & Epstein (2001), 

who required x very close to it/2 in order for their axisymmetric precession model to 

account for the observed precession of PSR 1828-11. (See also Rezania 2002.) Here, we 

~ 2 

also note that k 2 ~ e 2 f2_, where e 2 represents the deviation of the star from axisymmetry. 
Thus, nonaxisymmetric effects alone cannot introduce substantial harmonic structure in the 
phase residuals for small Cl_ either. 



The importance of spindown in the timing residuals is measure by the nondimensional 
parameter 

r ^ = -^h = ^hr - ^iooo^i (so) 

where P p = 1000P Pi iooo days is the precession period, Pq is the pulsar period in seconds, and 
t s d = lO 5 ^^ years. For PSR 1828-11, we have T s d ~ 844, so the pulsar spindown dominates 
the oscillatory terms. 



5.2. Pulse Arrival Times for \I B \ 3> C 

When \I B \ ^> C, the moment of inertia tensor is once again approximately 
axisymmetric, but neither 0_ nor + has to be small. Thus, the expansions in Appendices 
p.l| and D.2 are not applicable, and we shall have to solve the timing equation in a different 



way. For doing so, Eqs. (|(J2[ ) and flC^D prove to be useful. To keep the notation compact, 
we define the nondimensional parameter C = C/\I B \. 

Let us consider the two possible sign choices separately. For I B > 0, we use Eq. (|D24|) 
to evaluate b y . The problem is simplified since, from Eq. (|C2|), b + ~ C sin 26/2 <C 1. 
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Consequently, to first order in C, we find 

f ■ , Csin29 . . ,, 

o 3/ ~smxcos0 (sin r sm + cos x cos 0) , (ol) 

so pulses arrive when 

/ 1\ C sin 29 sin f 

0~ 2n+- 7T — ; 62 

V 2 / 2 sin % 

mapping from phase to time implies 

* dt L ( IN C sin 20 sin f C sin 2 # sin 2f , . 

-=— ~ ( 2ra + - ] 7r + . (63) 

7 V 2/ 2sinx 4cos% 

Taking account of pulsar spindown we find that the oscillatory part of the timing residuals 
is 



tt At osc ~ C 



sin 9 sin 2f sin 2$ sin f ( sin 6 1 cos 2f cos x sin 26 cos r 

4cosx 2sinx sd \ 16 2sinx 



• (64) 



Note that the phase residuals vanish as C — > even though the star still precesses. 
Moreover, there are oscillations at both uj p and 2u p whose amplitudes may be comparable, 
in agreement with observations of PSR 1828-11. 

For Ib < 0, we find, to order C, we find that pulses arrive when 

dt L ( 1\ (7sin26»cosT C sin 2 9 sin 2r 

— ~ 2n+- 7T — , (65) 

I \ 2 J 2 sin x 4 cos x 

and taking account of spindown results in oscilating timing residuals 



Q At osc ~ C 



sin 29 cos r sin 9 sin 2r ( sin 9 cos 2r sin 2^ cos x sin r 



+ r 



ml 



2sinx 4cosx V 16 2sinx 



(66) 



Once again, this involves oscillations at both uo p and 2uj p which can have comparable 
magnitudes. We see again that as C — >• 0, the oscillatory phase residuals disappear. 



5.3. Application to PSR 1828-11 

For PSR 1828-11, timing residuals appear to oscillate at both uo p and 2uo p with similar 
amplitude. The results of § |5.1| show that this situation is incompatible with small values of 
|f2_| if \Ib\ *C C. Moreover, we note that the results of § |5TT| continue to hold as \Ib\ — > 0, 
so we see that equal amplitudes at uj p and 2u p cannot arise from a model without magnetic 
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stresses, but with a triaxial crust, unless there is some fine-tuning of parameters (as in the 
axisymmetric model of Link & Epstein [2001], which requires x very close to vr/2). 



Thus, we focus on the strongly magnetic case, \Ib\ 3> C. In this case, Eqs. (J64j) and 
(|66|) show that it is possible for the phase residuals to oscillate with comparable amplitudes. 
The difference between the large and small \Ib\/C limits is that for small \Ib\/C, the 
amplitude of the observed timing residuals is determined solely by |fi_|, but at large \Ib\/C 
the amplitude is determined by C/\Ib\ primarily. Thus, by contrast to what we found for 
\Ib\ < C, small amplitude need not suppress the oscillations at 2ui p . 

For PSR 1828-11, we also know that T sc [ ~ 844 ^> 1, so let us approximate the 
oscillatory timing residuals further as 

CT s d sin 2 9 ( 16 cost \ , r . 

Q At osc ~ ^ cos2 r+— — (/b>0) 

tan U tan x J 




noAt „, c . | cos 2T _ (/B < o) . (67) 



The parameter 



16 

u = T——— (6? 



tan x tan ^ 

governs the relative strengths of the two harmonics in the timing residuals. If the precession 
is in the minimum energy state, 8 = x, then u = 1 for x — 76°. 

The timing residuals in Eq. flB?D vary between different minimum and maximum values. 
For Ib > 0, QoAt osc is minimum when cosf = —u/4, provided that |tanxtan#| < 1; for 
X = 6, this is so as long as x > 63°. Presuming this to be so, the minimum value is 

{iloAt osc ) min = — I 1 + — I . (69) 

The maximum is at cosf = 1, where we find 

(tt At osc ) max = CTadS * n (1 + u) . (70) 
16 

The ratio of maximum to minimum timing residual is 

1+u 



1 + u 2 /8 



(Ib > 0) ; (71) 



for u — 1, this ratio is 16/9 ~ 1.8, and the ratio is two for u — 2. Observationally, the 
timing residuals for PSR 1828-11 appear skewed toward positive values, with a maximum 
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about twice the magnitude of the minimum (see Fig. 1 in Stairs, Lyne & Shemar [2000]). 
Similarly, for < 0, the maximum value of f2 At osc occurs when sinr = — w/4, and the 
minimum occurs when sinr = 1; in this case, the ratio of minimum to maximum values is 

(il At osc ) max _l + u 2 /8 (J B < 0) • (72) 



l + u 

for u — 1, this ratio is 9/16 ~ 0.56, and the ratio is two for either u ~ 16.9 or u ~ —0.89. 
To the extent that we expect u > (manifestly so for \ = 8), and u ~ 1 (but not ^> 1) the 
observed residual arrival times appear to favor i# > 0, i.e. prolate magnetic distortions. 

Fig. ^] illustrates the residuals in the arrival times, At (dotted, left panel), period 
derivatives AP (solid, dotted) and period AP (solid, right panel) for 9 = x an d u = 5 for 
the prolate (upper) and oblate (lower) cases. The period an period derivatives are computed 
by differentiating Q At osc with respect to time: 

2tt d{Q At osc ) _ P 2 d{Q At osc ) 



AP(t) 



ttl dt 




27id 2 (tt At osc ) _ n _ ( P \ d 2 (Q At 
'o 

The agreement between these evaluations and the results plotted in Fig. 2 of Stairs, Lyne 
& Shemar (2000) is good, superficially, for the prolate model. Better agreement is seen for 
the variable period and period derivative than for the arrival times themselves; this may 
have been expected (e.g. Cordes 1993). 

Using the results plotted in Fig. [I], we can estimate the magnitude of C and the angles 
involved, even though these results do not constitute a true fit to the data, but just a 
plausible model. The curves for AP and AP resemble the observational results better, so 
let us focus on those. Fig. [I] was prepared for x — ® an d u = 5, which corresponds to 
X = = 60.8°. From Fig. 2 in Stairs, Lyne & Shemar (2000), we see that the maximum 
values of AP and AP are about 1 ns and 0.2 x 10~ 15 for PSR 1828-11. For the prolate 
model, which resembles the observations better, Fig. [I] implies a peak value of 460C ns for 
AP, and 2.QC x 10" 14 for AP. We therefore estimate G ~ 0.002 from AP, and C ~ 0.008 
from AP. Since we have not attempted true curve fits (i.e. by varying the parameters u, 9 
and C) we regard this as acceptable agreement, provisionally. 



6. Discussion 



Here, we have extended previous studies of precession of neutron stars to incorporate 
the effects of oblique magnetic fields. We have shown that if the magnetic stresses are large 
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1 2 3 1 2 3 

precession phase precession phase 



Fig. 1. — Results of evaluating the oscillating residual arrival time At and its first two 
derivatives AP and AP for u = 5 and x — ®- The top panels are for prolate models, the 
bottom for oblate. The left panels show At (dotted) and AP (solid), and the right panels 
shwo AP. The units of At are CT/Q , the units of AP are 2nCTu p /Q 2 = CTP 2 /P P 
(approximately l.QOCfis for PSR 1828-11), and the units of AP are 27rCTwJ/fig = 
2ttCTP 2 /P p 2 (approximately 1.17(7 x 1(T 13 for PSR 1828-11). 
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enough, then steady rotation is unlikely, and the neutron star must precess. Moreover, 
even when the magnetic stresses are relatively weak, so that steady rotation is possible 
irrespective of the obliquity of the magnetic field, the minimum energy state is not the 
steady state. Thus, even in this case, the neutron star will precess. We argued, in §[|, 
that the minimum energy, precessing state is a local energy minimum that applies at fixed 
angle between the magnetic and rotational axes. On a longer timescale, we would expect 
the star to seek its global energy minimum, which should correspond to either aligned 
or perpendicular magnetic and rotation axes, and no precession. We might expect short 
timescale dissipative effects to drive the system toward its local minimum, and that the 
global minimum is only achieved on somewhat longer timescale, perhaps as a result of 
electromagnetic spindown torques (e.g. Goldreich 1970). 

The effective moment of inertia tensor of a neutron star with an inclined magnetic field 
is inherently triaxial. Consequently, the precession is periodic but not sinusoidal in time. In 
general, the solution for the rotational angular velocity of the star can be expanded in a 
Fourier series involving harmonics of the precession frequency. We have shown that at least 
the first few terms in such an expansion can have comparable magnitudes provided that the 
interior magnetic stresses are not very small. 

The condition that magnetic stresses play an important role is that the magnetic- 
induced distortions are comparable to or larger than the distortions of the stellar crust. 
For precession periods of order years, the implied magnetic stresses exceed those expected 
from the classical Maxwell stress tensor, evaluated using the inferred dipole magnetic 
field strength, by a couple of orders of magnitude. However, if the interior of a neutron 
star contains a Type II superconductor, or else is a normal conductor but possesses large 
toroidal magnetic fields, the magnetic stresses are larger, and the implied distortions can be 
of the right order of magnitude (Jones 1975, Easson & Pethick 1977, Cutler 2002). Thus, 
the observation of neutron star precession can be taken as indirect evidence for enhanced 
magnetic stresses, due to either Type II superconductivity, or large toroidal fields. 

We postpone detailed application of the ideas set forth here to PSR 1828-11 to another 



paper (Akgun, Epstein & Wasserman 2002). However, in §:k3 we argued that only a model 
with \Ib\ ^> C can lead to time residuals that oscillate with comparable amplitude at 
both uj p and 2u p . In this case, the amplitude of the observed time residuals is set by the 
dimensionless ratio C = C/\Ib\, not by the tilt of the angular velocity vector away from 
any principal axis, which may be large. 

By contrast, if the stellar distortions associated with magnetic stresses are smaller than 
those associated with the crust, then Eqs. ( |57D and (|58|) show that the timing residuals 
oscillate predominantly at uo p . Oscillations at 2uo p would be down by factors ~ \&-\, where 
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fl is given by Eq. (p9[), and is ~ \Ib\/C in the minimum energy configuration. Moreover, 
we argued, in § pTT| and § |5^ , that precession of a triaxial crust alone would probably not, at 
small precession amplitude, be capable of producing oscillations of comparable magnitude 
at both Up and 2u p , because the precession amplitude is proportional to |f2_| and the 
oscillations at harmonics of u p are suppressed by factors ~ \&-\- Thus, a solution in which 
the crustal distortion is responsible for precession is unlikely to explain the data on PSR 
1828-11. The model in which magnetic stresses dominates is not merely precession of a 
triaxial body because small amplitude phase residuals can arise even when |f2_| is not 
small. 

For a precessing star in which the distortion is due to magnetic stresses primarily, the 
precession period is 

P I 492 days 

P ~ ~ (741 

p 3J B cosx P cos x{BH) 27 / '5 ' 1 ; 

where the value given is for PSR 1828-11. This is similar to the observed period, about 1000 
days, provided that (3 cos \ is n °t very small. Thus, if % ~ 60° the timing residuals could 
oscillate with about the right period, amplitude and relative importance of the fundamental 
precession period and its first harmonic. In this case, we note that the condition that there 
is no steady state, | sin2x| > C can be satisfied: for \ — 60°, for example, sin 2% = V3/2 
whereas we estimated that C ~ 0.001 — 0.01 in § |5.3| . This possible explanation of the 
timing residuals for PSR 1828-11 only works if the magnetic stresses in this star are ~ 200 
times larger than would be indicated by its dipole magnetic field. Thus, we conclude that 
either the interior is a Type II superconductor, or is a normal conductor with a toroidal 
field whose strength is ~ 10 14 G. Otherwise, the expected magnetic stresses are far smaller 
than is needed for this solution to apply. 



We also noted, in § |5.3| , that the ratio of the minimum and maximum timing residuals, 
and the shape of the variation of AP and AP seen in PSR 1828-11 appear to favor a model 
in which Is > 0, so that magnetic distortions are prolate. Prolate distortions would arise 
naturally from the stresses due to a toroidal field, with or without Type II superconductivity 
(Cutler 2002), but may also result if magnetic flux tubes have been transported outward in 
the core and accumulate at its outer boundary (Ruderman, Zhu & Chen 1998, Ruderman 
& Chen 1999), which could "pinch" the interior. 

Crustal distortions are still needed in order for the precessional amplitude to be 
nonzero. In fact, to be more precise, the precessional amplitude depends on the component 
of the moment of inertia tensor of the crust that is not symmetric about the magnetic axis. 



This can be seen directly from Eqs. (^4|) and (66), which show that the oscillating time 



residuals vanish as C sm9 — > 0, where 6 is the angle between the magnetic axis b and, in 
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this axisymmetric distortion model, the symmetry axis of the relevant crustal deformation, 
k. Thus, although all neutron stars precess as a consequence of their magnetic stresses for 
| is | 3> C in the picture advanced here, only those with sufficiently large nonaligned crustal 
deformations would have discernible oscillations of their timing residuals. In this sense, 
PSR 1828-11 may be special. 

Although we have treated the basic physics of precession of an oblique rotator in some 
detail, we have not treated several effects that might play significant roles. We have not 
explicitly included either vortex line pinning or vortex drag. Link & Cutler (2002) have 
argued that vortex lines can unpin globally at large enough precession amplitude. For 
\Ib\ C, as is required to explain the timing residuals in PSR 1828-11, the precession 
amplitude is ~ sinx or cosx (depending on the sign of Is), which is not small, so global 
unpinning is expected. Thus, we may expect that vortex lines are unpinned in neutron stars 
with magnetic fields that are strong enough to have \Ib\ 3> C except for \ very close to 
either zero or tt/2, depending on the sign of Jg. For these, vortex line drag, if weak enough, 
may simply serve to bring the neutron star superfluid into corotation with its crust, and 
drive the rotating star toward its minimum energy state. For large \Ig\/C, we have seen 
that precession is required, so weak vortex drag or other forms of dissipation need not 
prevent precession. For small values of \Ib\/C, the precession amplitude is given by Eq. ( |59"D 
and can be very small; in the minimum energy state, Eq. (|59| ) implies |Q_| ~ |I_Bsin20|/C. 
In this case, it is possible that precession cannot overcome pinning forces, as discussed in 
Link & Cutler (2002), and so precession does not occur. 

Theories of pulsar glitches involve the pinning, unpinning and repinning of crustal 
superfluid vortex lines (e.g. Anderson and Itoh 1975; Alpar et al. 1981, 1984a,b, 1993; Link, 
Epstein and Baym 1993). As we have seen, for large \Ib\/C, precession amplitudes are 
large, and vortices may be expected to unpin, but it is possible for vortex lines to remain 
pinned in neutron stars with C 3> \Ib\- Thus, there could be a dichotomy between pulsars 
that glitch {C ^> \Ib\) and those that precess (C <C |-Tb|)- If, in the course of a glitch, all 
crustal superfluid vortices were to unpin, then the star might precess briefly. Perhaps that 
explains the detection of damped, quasisinusoidal timing residuals in the Vela pulsar after 
and perhaps before its Christmas 1988 glitch (McCulloch et al. 1990). 

We have kept the problem of precession of an oblique rotator as simple as possible by 
considering what happens when the magnetic field is axisymmetric about some axis, and 
the crustal distortions are also axisymmetric, but about a different axis. More realistically, 
both of these simplifying assumptions are likely to be violated. Most likely, the crust is 
not axisymmetric. When magnetic stresses dominate, we do not expect including intrinsic 
crustal asymmetry to alter the results found here qualitatively, since the effective moment 
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of inertia is already triaxial here. Triaxiality of the crust, in the limit of rather small crustal 
distortions, would simply rotate the principal axes slightly. Furthermore, the magnetic field 
may have a more complicated structure than we have assumed. A substantial quadrupolar 
component would presumably render the contribution to the inertia tensor from magnetic 
stresses alone triaxial. We shall consider these complications elsewhere. 

Although we have included the spindown torque in our evaluations of timing residuals, 
we did not include near zone electromagnetic torques (Good & Ng 1985, Melatos 1997, 1999, 
2000). The principal effect of such torques would be to renormalize the moment of inertia 
tensor of the star. Near zone torques can play a role similar to the magnetic distortions 
considered here, but are smaller by a factor ~ (H/B)(Rc 2 /GM), which is non-negligible 
even if H = B. However, we note here that the large magnetic distortions we propose would 
presumably apply to the spindown of magnetars and anomalous X-ray pulsars, in much the 
same fashion as proposed by Melatos (1999, 2000). We shall pursue this idea elsewhere. 

We have also ignored motions of the fluid and crust of the star apart from rigid 
rotation. Mestel and collaborators (Mestel & Takhar 1972, Mestel et al. 1981, Nittman & 
Wood 1981) have pointed out that the equilibrium in a fluid star with oblique magnetic 
field must involve fluid motions with velocities ~ (f2 2 i? 3 /GM)uj p R. These distort the 
stellar magnetic field by a fractional amount ~ Q 2 R 3 /GM over the precession period 2n/ujp 
(Mestel & Takhar 1972, Mestel et al. 1981, Nittman & Wood 1981). Similar motions might 
arise in the crust, but with magnitudes ~ (Q 2 R 2 /c 2 )u p R, where c t is the sound speed for 
transverse waves. The expected amplitude of the resulting magnetic wander is ~ Q 2 R 2 /c 2 . 
Although slow, these displacements cause the magnetic field of the star to oscillate about 
its undisturbed, axisymmetric state, and might influence the long-term behavior of the 
observed spindown. We shall investigate whether there is any long-term observational 
signature of these motions elsewhere. In addition, it is likely that magnetic and rotational 
deformations of the core must also deform the crust, as they exert pressure on its inner 
boundary. We would expect the magnetic deformation of the core to promote crustal 
deformation symmetric about b, which would not lead to observable precession, but the 
rotation-induced deformation need not be symmetric about b, and should be substantial. 
An important question left unanswered here is whether the crustal distortion, dj, that 
leads to detectable precession is relatively steady, or is simply due to a seismic fluctuation. 
We shall explore the important issue of crustal deformations elsewhere. 

Finally, we emphasize that any neutron star with strong enough core magnetic stresses 
ought to precess, but we may not be able to detect their precession because their magnetic 
axes can still rotate more or less uniformly. This is because, at small values of C/\Ib\, 
the neutron star precesses almost exactly about its magnetic axis, which therefore rotates 
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almost uniformly as seen in the inertial frame. Although the precession may not be 
detectable readily from timing residuals for most pulsars, gravitational radiation amplitudes 
would be larger than would arise without enhanced internal magnetic stresses (e.g. Cutler 
Sz Thorne 2002, Cutler 2002). The distortions required for PSR 1828-11 are still smaller 
than would be needed for detection by LIGO, even if it were spinning faster (Brady et al. 
1998). If there are young, highly magnetized neutron stars rotating rapidly, they would be 
the brightest emitters of gravitational radiation. Such objects have been hypothesized to be 
the sources of the highest energy cosmic rays (Blasi, Epstein & Olinto 2000, Arons 2002). 

Partial support for this work was provided by a grant from IGPP at LANL. I thank T. 
Akgun, J. Cordes, R. Epstein and B. Link for comments. 



A. Inequalities Among Eigenvalues 



For I B > 0, we can rewrite Eq. ( |5T] ) in the form 
3 



6L 



from which it follows that 



C-I B ± J(C - I B ) 2 + 2C/ B (1 - cos 26) 



(Al) 



51+ > -(C-I B 
51- < \{C-I B 



\C-Ib\) 
\C-Ib\) ■ 



(A2) 



Eq. implies that 61+ > 3(C - I B ) > and 5J_ < if I B < C, and 5I + > and 

51- < 3(C - I B ) < if I B > C. Thus, if I B > 0, then 5I + > and 51- < irrespective of 
whether I B > C or I B < C. Thus, for I B > 0, i+ > J 2 > I— 



(A3) 



For Ib < 0, we can rewrite Eq. (|3TD in the form 
3 



5U 



C+ \I B \ ± J(C- \Ib\) 2 + 2C\I b \(1 + cos26) 



from which it follows that 



SI+ > -(C+\I B \ + \C-\I B \ 



51- > ~(C 
~ 2 K 



l B\ 



\c 



l B 



. 



(A4) 



Eq. ( |A4j ) implies that 5I + > 3C if \I B \ < C and 51 + > S\I B \ if \I B \ > C. Thus, we see that 
for I B < 0, J + > I > J 2 . 



B. Approximate Results for \I B \ <C C 



For \I B \ <C C, we can expand these results to find 

1 B cos 2 7|(5 cos 4 - 3 cos 2 1 



sin a + 



COSCTj 



Using Eq. 



-b 

-b, 
51, 
61- 



sin 

cos 9 



1 + 



C 



2C 2 



1 J B sin 2 9 I B (5 sin 4 0-3 sin 2 , 



C 



2C 2 



3C 



Ml + cos 20) i| sin 2 20 
1 - — — — + , —„ h ■ ■ 



2C AC 2 
3/^(1 -cos 20) 3/|sin 2 20 



4C 2 



+ ■ 



we also find that 
A 



-(o) f2_ issm20 
" + + 2C + ' 
^(o) _ Of J B sin 29 



3£/ 
2/ I 



3C 



2C 
I B cos29 



C 



+ 



2/ 



C 



(0) 



where + = cos(x — 9) and 0_ = sin(x 



■ + 
+ • 



C. Approximate Results for \I B \ C 



For \I B \ ^> C, we find 

7/ 



1 ± 



61+ 
61- 



~ 1 ± s B 

3|Jb| 
2 

3|/b| 



, C 2 sin 2 29 

2?| + 
C(l -sb cos 20) C 2 sin 2 20 

C(l + s B cos29) C 2 sin 2 20 
-1 - s B + — rrn — + 



2/ 2 



where s B = I B /\I B \ is the sign of i^. For I B > these results imply 



sin ct-j 



C 2 sin 2 20 
-b ~ 1 h 
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C sin 26 

cosa+ = 6+ ~ — — + ••• 

~ C cos x sin 26* 
fi+ ~ sin x + ^7 + • • 

A C sin x sin 26* 
S2_ ~ -cosxH — h 

3J B / Ccos20 \ 

A ~ -=-[1 ; + ■■■] 

21 \ I B ) 

_ 2 A + A Csin 2 6> 



e ~ 



A - A / 



B 



2AL cos x , no , 
u p ~ = , (C2) 



and for Ib < the same results imply 

C sin 20 

sm<7 + = — b ~ — — — — : h • • • 

2j Jb| 

C 2 sin 2 2# 
cos a + = o + ~ 1 -7- 1 

- C sin x sin 26* 

n+ ~ cos x + — -sm — + 

- . C cos y sin 26* 
ft„ ~ smx ;r— + 

z \ 1 b\ 

3\I b \ ( Ccos29 
A ~ -U^- 1 



2/ V J 



2 A -A Csin 2 fl 

e ~ - — — — : — : h • • • nonumoer (Co) 

2A |/ B | 

2AL V / e 2 + ri+ 2ALcosx 
cj p ~ Sf— i ~ = — - . (C4) 



D. Timing Solution 

To find pulse arrival times, we need to determine the motion of b in the inertial 
reference frame of the observer. To do this, we need the Euler angle rotation from the 
rotating frame of reference to the inertial frame; for an arbitrary vector V this is (see e.g. 
Goldstein Eq. [4-47]) 

V x = V j (cos ip cos — cos a sin -0 sin 0) — V " 2 (sin -0 cos + cos a cos ip sin 0) + V§ sin a sin 
V y = V j (cos -0 sin + cos ct sin -0 cos 0) — V 2 ( s - n "0 sin — cos a cos -0 cos 0) — V§ sin a cos 
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V z = Vj sin a sin ^ + sin a cos ^ + V§ cos a , (Dl) 

where a, 0, ^ are the Euler angles defined in Fig. 47 of Landau & Lifshitz, §35, except that, 
to avoid confusion with our definition of 9 as the angle between k and b, we label their 
Euler angle 9 as a. We assume that L is along the e z direction in the inertial frame. We 
can then determine the two angles a and i/j from the equations 

Lg = L cos a 

Lj = L sin a smip 

Lrj = L sin a cos ip ; (D2) 
the third Euler angle is not determined by these relations, but can be found from 

dt - {iM + i?ni) ■ (D3) 

The choice of axes ej § § i n the rotating frame of reference is somewhat arbitrary, and we 
shall make three different choices below, as the situation demands. 



D.l. Pulse Arrival Times for I B > 

D.l.l. I B >0 and k 2 < 1 

For I B > and k 2 < 1, we choose Vj = V_, = V 2, and V 5 = V + ; then Eqs. (^) 
imply (Landau & Lifshitz Eq. [37.15]) 



cos a 



£l+(l + A) dn(r) 
1 + 2A(2f2+ - 1) -\ 
cn(r) 



A 2 



Cl + dn(r) 



(A + A )(l-A) cn(r) (A + A c 



sn(r) \\ 2A(1 - A ) sn(r) 



2A 



(D4) 



where the approximations are to leading order in A. The remaining Euler angle <fi evolves 
according to (Landau & Lifshitz, Eq. [37.16]) 



d(f) 
~dt 



L 
1 
L 
1 



A + A + (A- A ) sn 2 (r) 



1 - A)(A + A ) + (A - A )(l + A) sn 2 (r) 
A[A + A -(A-A ) sn 2 (r)]' 



1 + 



A + A + (A - A ) sn 2 (r) 



(D5) 
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which is uniform up to corrections ~ A. Note that in perfect axisymmetry, fl z and d<p/dt 
are independent of time, another difference between the magnetic case, which is inherently 
triaxial, and precession with an axisymmetric crust. When the star is nearly axisymmetric, 

^£[l + A-(A-A )sn 2 (r)] ; (D6) 

thus, d(f)/dt oscillates at twice the precession frequency in this limit. Associated with 
the time development of dcp/dt would also be variability of the pulsar spindown, at even 
harmonics of the precession frequency. The amplitude of the main variation would be of 
order (A — A )P /P P ~ e(P /Pp) 2 , where e ~ (A — A )/2A is the fractional deviation from 
axisymmetry; successive harmonics would be smaller by powers of e. For PSR 1828-11, 
we would have (A — Ao)Pq/P p ~ 3 x 10~ 16 e. For comparison, electromagnetic spindown 
produces oscillations with an amplitude ~ |Q_|Po/ r sd ~ 1-3 x 10~ 13 |f2_|, where T s d is the 
spindown timescale, which is considerably larger (Link & Epstein 2001). 



Applying Eqs. ( Pl|) to b = b + e + + 6_e_, and using Eqs. 



1 + 



sn2(r)(A-A )(l+A) 



sn t smi 



(A+A )(l-A) 

6 + n_(l - A)cos0 
l + 2A(2ft+ - 1) + A : 



2A(1 - A r 



we find 



+ A) dn(r) cn(r) cos . 



^(A + A )(l-A) ^/i + 2 A(2n+-l) + A= 



1 + 



sn 2 (r)(A- A )(l + A) 
(A + A )(l-A) 



(D7) 



No approximations have been made in Eq. (p7|) ; in fact, there is also no explicit reference 
to magnetic distortions here, and so these results apply to triaxial stars in general. We note 
here that 

1/2 

If 'ens 'It) - In) 



1 

71 



C 2 - 2I B C cos29 + I 2 



B 



b = = 



1 

V2 



(Ccos20-I B ) 



C 2 - 21 B C cos 29 + I 2 B 



1/2 



(D8) 



and in the minimum energy state we can take 9 ~ \ to lowest order in distortions. When 
C is large compared with I B , these reduce to b + ~ cosx and 6 ~ — sin^. 

Pulse arrival times are found from the condition that b y = if we assume that the 
observer is in the x — z plane. Let us simplify the notation by introducing the parameters 



e 2 = 



(A-A )(l + A) 
(A + A )(l-A) 
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q± 



1 ± A 

1 + 2A(2f2+ - 1) + A 2 ' 



(D9) 



e 2 measures the departure from axisymmetry, and q± are one to lowest order in the 
distortions. We can simplify further by defining S± = q±£l±; then it follows that 

-2 -2 

S + + £_ = 1. In terms of these, the pulse arrival times are the solution to 



b_ [vl + e 2 sn(r) sin + S + dn(r)cn(r) cos 0] 
y/l + e 2 sn 2 (r) 

we can rewrite this as 

- r) = [sinr — y/l + e 2 sn(r)] sin0 



b + X cos 0\/ 1 + e 2 sn 2 (r 



(D10) 



cos 



+ I cost - S + cn(r)dn(r) + + e sn ( r )] ? COSl 



(Dll) 



Define the functions 



e 2 F s {r) 
e 2 FMr) 



sinr — vl-j- e 2 sn(r) 



"CV -7 J = cost — cn(r)dn(r) ; 
then the pulse arrival times are the solutions of the equation 



(D12) 



cos 



r) = e 2 F s (T) sin0 + e 2 F c (r)S + + (1 - S+) cost 



~XJ_[1 + e 2 sn 2 (r)] \cos> 
6_ 



(D13) 



In general, we can solve Eq. ( p!3| ) numerically, but a good analytic approximation can be 
found if the RHS is small. 

The terms on the RHS of Eq. ( pi3| ) will be small if the deviation from axisymmetry 
is small and the rotation is nearly aligned with one of the principal axes of the inertia 
tensor. If we assume that RHS is small, then we can develop an approximate solution by 
first substituting = — r + r. If we let r\ = — r, then we find the equation 



cos?7 = sin 77 
where 



e 2 Fs(r) cos r — F(r) sin r + cos rj e 2 Fs(r) sin r + F(t) cos r 



(D14) 



F(t) = e 2 F c (r)E + + (!-£+) cost + ^£_[l + e 2 sn 2 (r)] . 



(D15) 
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We see that F(r) is small only if both e 2 «l and 1 — S + ~ 1 — f2 + <C 1, so the expansion 
implied here involves a pair of small parameters. To zeroth order the solution to Eq. ( pi4 ) 
is rj = (2n + 1/2) it; writing the solution in general as rj = (2n + 1/2) 7r + S we find 



tan 5 



[e 2 F s {r) cost — F(t) sinr] 
— [e 2 Fs(r) sin r + F(r) cos r] 



or, for small values of 8, 



8 ~ [F(t) sin t - e 2 F s (r) cos r)] 1 + e 2 F 5 (r) sin r + F(r) cos r + 



(D16) 



(D17) 



It should be apparent from Eq. (p!7|) that arrival times should exhibit oscillations of 
decreasing magnitude with frequencies ru p , where r is an integer. Rigorously, we should 
expand the functions F(t) and Fs(r) in powers of e 2 as well, and the hiden functions dn(r), 

A 2 

cn(r) and sn(r) should be expanded in k 2 (which is small if 1 — Q + is small: see Eq. [[4!| ). 

For the triaxial case we have to consider, in addition to 5, the mapping between t and 
0, which is complicated by the fact that d<p/dt is time dependent. Eq. ( |D5| ) implies that, to 
lowest nontrivial order in e 2 , 



--j- ~ 1 + A(l - e 2 ) + e 2 A cos2r 



which integrates to 



and therefore 



1 + A(l - e 2 



* dt L e 2 sin 2r 

7 4in+i 



1 - A(2|n + | - 1 + e 2 



* dt e 2 sin 2r 
? 4IO.I 



(D18) 



(D19) 



(D20) 



The integral includes the effects of sinusoidal variations in pulsar spindown, which may be 
evaluated using the formulae in § fOI. Thus, the final solution for arrival times is 



1 - A(2|n+| - 1 + e 2 



t dt L 



2n + - ) 7T + 5 - 



e 2 sin 2r 
4IOJ 



(D21) 



We see that the variation of d<fi/dt with time may introduce variation in pulse arrival times 
at twice the precession frequency, a possibility that only arises when the moment of inertia 
tensor is nonaxisymmetric. Note that in the expansion in Eq. ( piS| ) we only kept terms up 
to ~ e 2 . Higher order terms would introduce further variability (presumably at 4uj p , 6uj p 
and so on). 
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D.1.2. I R >0 and k 2 > 1 



For k 2 > 1, we choose Vj = V + , = V 2 and V§ = V_. We then find that 



cos a = (1 — A)dn(f) 



1 - O 



\ 1 + 2A(2f2 + - 1) + A 2 



t&nip = sign(0. 



cnr 



(1 + A)(A- A ) 



dt 



L 
1 



sn(f)^ 2A(l-Ao) 
A - A + (A + A )sn 2 (f) 



(1 + A) (A - A ) + (1 - A) (A + A )sn 2 (f ) 



(D22) 



For this case, we define S± as before, and 



g 2 



(A + A )(l-A) 
(A - A )(l + A) 

in terms of these variables, we find that 

- 6 + [sign(r2_)v / l + e 2 sn(f) sin0 + E_dn(f) cos0] 

1 + e 2 sn 2 (f) 



(D23) 



b_£ + cos<Wl + e 2 sn 2 (r) , (D24) 



where b± are given by Eq. ( |D 8| ) . Pulse arrival times are found by solving b y = 0. The 
equations are exactly the same as for k 2 < 1 except for the replacements e 2 — > e 2 = 1/e 2 , 
b — > sign(f2_)6 + , b + — > b . S + — > and S_ — > E + . Thus, if we define rj — <fi — f, and 
also 



e 2 Fs(f) 



sin r — vT+ e 2 snl r ) 



e Fc(t) = cost — cn(f)dn(f) 



F(f) = e 2 F c (f)|S_| + (l- |S_|)cosr + sign(n_)^S + [l + g 2 sn 2 (f)] , (D25) 



then we find 



cos?7 = sin r/ 



e Fs(t) cos f — .F(t) sin r +cosr/ e Fs(f) sinf + F(f) cos-f . (D26) 



As before, the approximate solution is r\ = (2n + l/2)7r + 5, where 

[e 2 Fs(f) cos f — F(f) sinf] 



tan 5 



1 - [g 2 F 5 (f ) + F(f ) cost] 



(D27) 
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which becomes 



5 ~ [F(t) sin r — e F s (r) cos r] [1 + e F s (t) sin r + F(r) cos r + 



(D28) 



for small 5. This approximate solution is valid as long as |f2_| ~ 1 and e 2 <C 1. The 
restriction to small |f2_| may seem a bit strange, until we recall that for 1b > and 
k 2 = 1/k 2 < 1, the spinning star is somewhat prolate rather than oblate. Thus, simple 
solutions are expected for small amplitude around the axis of smallest moment of inertia in 
this case. 

The mapping from pulse phase to pulse arrival times requires the connection between 
7] = cf) — f and t. From the last of Eqs. (|D22|) we find that 



--j- ~ 1 - A(l - e 2 ) - e 2 A cos 2f , 

Li dt 



which integrates to 



Afl 



dt L e 2 sin 2r 



4\ i-n 



~~T ' 



(D29) 
(D30) 



when we use cu p ~ 2ALy 1 — Sl + /L Thus we find that 
?7 = 0-f~ 1- A ^1 + 2^/l - f2+ - e 2 ^ 

so that pulses arrive at times 

1 - A ( 1 + 2a/ 1 - &1 - e 2 I / -V 1 ^ ( 2/' 



* dt L e 2 sin 2f 



4a/1 - ft 



1 \ r e 2 sin 2r 



(D31) 



(D32) 



As before, we see that to lowest order in e 2 , the mapping from phase to time introduces 
variation at 2u p . Retaining higher order terms would introduce further variability at 4u; p , 
etc. 



D.2. Pulse Arrival Times for Ib < 

For I B < 0, Vj = V 2 and V§ = V_, and Eq. ( p4|) is replaced by 

dn(r) 



cos a 



tan^ 



[A + A 2 + A(l + A )(2ft + - !)](! + A) 
\ (A + A ) [1 + 2A(2S7+ - 1) + A 2 ] 



cn r 



2Afl - A r 



sn(r)\ (1-A)(A + A C 



(D33) 
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and Eq. (|D5|) is replaced by 

d<p L 
~dt ~ 7 

In this case, we define 



2A + (A - A)sn 2 (r) 



2A(1 - A ) + (A - A)(l + A)sn 2 (r) 



(D34) 



(1 + A)[A + A 2 + A(l + A )(2O + - 1)] 
\\ (A + A ) [1 + 2A(20+ - 1) + A 2 ] 



2A(1- A )(l- A) 



\ (A + A )[l + 2A(2f2 + -l) + A 2 ] 
(A -A)(l + A) 
2A(1 - A ) ' 



(D35) 



then 



6_[cn(r) sin0 — S + a/1 + e 2 dn(r)sn(r) cos< 
1 + e 2 sn 2 (r) 



- b+S_ cos cWl + e 2 sn 2 (r) , (D36) 



where b± are given by Eq. (|D8| ). Pulse arrival times are determined from the condition 



b y = 0: 



sin ffl — t) 



[cos r — cn(r)] sin <p 

+ js+Vl + e 2 sn(r)dn(r) - sinr + ^^-l 1 + e2sn2 ( r )]| cos< 



(D37) 



Define the functions 



e 2 Gcij) = cosr — cn(r 



e 2 G, 



vl + e 2 sn(r)dn(r) — sinr 



G(t) = S + e 2 G 5 (r)-(l-S + )sinr+^S_[l + e 2 sn 2 (r)] 



(D38) 



then the pulse arrival times are solutions of the equation 

sin r\ = [e 2 Gc(' r ) cos T — G(t) sin r] sin rj + [e 2 Gc(r) sin r + G(r) cos r] cos rj , (D39) 
where rj = <fi — r. The solution is rj — 2im + 5, where 

e 2 Gc{T) sin r + G(r) cos r 



tan 5 



1 — [e 2 Gc(j) cos r — G{r) sin r 



(D40) 
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or, for small 5, 



S ~ \e 2 G c {r) sin r + G(r) cos r] [1 + e 2 G c (r) cos r - G(r) sin r] . 



(D41) 



Once again, it is necesary that both e 2 C 1 and J 1 — f2 + <C 1 for the approximation to be 
valid. 

For this case, the mapping from phase to arrival times that results from expanding Eq. 
( p34| ) leads to exactly the same result as for Ib > and k 2 < 1, namely 



l-A(2f2 + -l + e 2 )]^ 



«■* dt L e 2 sin 2r 
-=- + : * : 



(D42) 



In this case, then, pulses arrive when 



1 - A (2ft+ - 1 + e : 



* L 
o ~T~ 



2n7r + 5 



e 2 sin 2r 
4|A7i 



(D43) 
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